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Extending the single photon transmission distance is a basic reqnirement for the implementation 
of quantum communication on a global scale. In this work we report the single photon exchange 
from a medium Earth orbit satellite (MEO) at more than 7000 km of slant distance to the ground 
station at the Matera Laser Ranging Observatory. The single photon transmitter was realized by 
exploiting the corner cube retro-reflectors mounted on the LAGEOS-2 satellite. Long duration of 
data collection is possible with such altitude, up to 43 minutes in a single passage. The mean 
number of photons per pulse (Usat) has been limited to 1 for 200 seconds, resulting in an average 
detection rate of 3.0 countsis and a signal to noise ratio of 1.5. The feasibility of single photon 
exchange from MEO satellites paves the way to tests of Quantum Mechanics in moving frames and 
to global Quantum Information. 


Introduction - Quantum Gommunications (QC) are 
necessary for tests on the foundation of Quantum 
Physics, such as Bell’s inequalities violation m , entan¬ 
glement swapping Hig and distribution [B], and quan¬ 
tum teleportation [7] . Moreover, the transmission of light 
quanta over long distances is crucial for the realization of 
several Quantum Information protocols, as quantum key 
distribution (QKD) [SUIT)], quantum authentication [11] 
and quantum digital signature [T^ . 

One of the main challenges of QC consists in extend¬ 
ing the length of the communication channel, which on 
ground is nowadays limited to few hundreds of kilome¬ 
ters, both for free space |7] and optical fibers [13]. To 
overcome these limitations, in the last decade the feasi¬ 
bility of a quantum channel from an orbiting terminal 
has been explored both theoretically [IMIS and exper¬ 
imentally by exploiting the Ajisai m and Champ |21] 
satellites. Notably, recent transmission of qubits beyond 
Earth atmosphere [52] proved the feasibility of QC from 
several low Earth orbit (LEO) satellites, within 1500 km 
of altitude with respect to Earth surface and with typi¬ 
cal passage duration of a few minutes. However, propos¬ 
als for fundamental tests of Quantum Mechanics with 
moving frames, involving Special and General Relativ¬ 
ity are based on link distances exceeding the LEO scale 
[551124] . Moreover, satellites with orbits higher than LEO 
are interesting for the implementation of unconditionally- 
secure information protocols for future global navigation 
satellite system (GNSS) constellations [26] or for the re¬ 
alization of permanent links with geostationary (GEO) 
satellites. 

In this perspective, this work aims at the expansion 
of QC with higher satellites by realizing the first single 
photon exchange spanning a distance exceeding 7000 km, 
with a mean photon number per pulse leaving the satel¬ 
lite (psat) lower than 1 for an extended part of the link 
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FIG. 1. (Color online) Top: schematic of the two-way pro¬ 
tocol. A pulsed laser beam is directed toward a satellite 
equipped with OCRs and redirected toward the ground sta¬ 
tion. The two satellite trajectories are drawn to reproduce 
the distances in scale of Ajisai and LAGEOS-2 from the Earth 
surface. Bottom: measured radial velocity vr in natural units 
P = vr/c. The corresponding Doppler effect on the pulse sep¬ 
aration is schematized on the three boxes: when /3 < 0 (/3 > 0) 
the pulse separations at the receiver are smaller (larger) than 
the transmitted pulse separations. The Doppler effect affects 
in the same way the 10 Hz SLR beam and the 100 MHz 
beam used for single photon generation. 


and with the single photon return signal above the noise. 
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Moreover we will show that, with an upgrade of the de¬ 
tector on ground with commercially available parts, it 
would be possible to achieve signal to noise ratio (SNR) 
suitable for quantum communication up to 23000 fcm, 
corresponding to the distance to GNSS. 

Experimental realization - As sketched in the top panel 
of Fig. to simulate a single photon source in orbit 
we exploited a train of light pulses directed towards the 
LAGEOS-2 satellite and retroreflected by the satellite 
corner cube retroreflectors (GGRs). Such pulses are at¬ 
tenuated in the uplink by the combined effect of telescope 
divergence, optical diffraction, atmospheric absorption 
and atmospheric turbulence |25j . The mean photon num¬ 
ber /Xsat of the retroreflected pulses has been estimated 
by the link budget equation. Because of OCRs backre- 
flection capacity, the ground station is automatically in 
the illuminated cone on ground. 

The pulses are directed towards the satellite by a track¬ 
ing telescope that is also used as single-photon receiver. 
The transmitter and receiver setups share the same tele¬ 
scope Coude path and are combined by a beam splitter 
(BS). The train of upgoing pulses is generated by a mode¬ 
locking Nd:YV04 laser oscillator, operating at 1064 nm 
and with repetition rate stabilized at 100 MHz by an 
atomic clock. By means of a PPLN nonlinear crystal, 
we obtained pulses at the wavelength of 532 nm, with 
1.1 nJ of energy and a duration of 100 ps FWHM. The 
beam size and divergence are corrected with a Galilean 
beam expander before entering the Coude path of the 
telescope. Finally, the pulses are directed towards the 
satellite by the 1.5 m diameter primary mirror of the 
Matera Laser Ranging Observatory (MLRO) of the Ital¬ 
ian Space Agency, located in Matera, Italy. 

The optical setup for the detection of the single pho¬ 
tons from the satellite is composed by a focusing optics, 
field-of-view (FOV) control, filtering and a single photon 
detector, based on a photomultiplier tube (PMT), Hama¬ 
matsu H7360-02. The PMT TTL output is time-tagged 
with a 81 ps time resolution by a fast time-to-digital con¬ 
verter (TDC), QuTAU from QuTOOLS. Two shutters 
are used to alternate the transmission and reception pe¬ 
riods and to avoid the exposure of the PMT to stray 
upward light. The FWHM of the receiver FOV has been 
set to 100 prad for reducing the gathering of background. 
In addition, a spectral filter of 3 nm FWHM transmis¬ 
sion band around 532 nm has been used to reduced the 
stray light counts below the 50 counts/s intrinsic dark 
counts of the PMT. 

Temporal synchronization - To discriminate photons 
received from the satellite from the uncorrelated back¬ 
ground, detected events are temporally filtered accord¬ 
ing to the instantaneous round trip time of light. Photon 
expected-times-of-arrival Eef at the receiver are not pe¬ 
riodic during the satellite passage, as they depend on the 
actual distance from the ground station and on the ra¬ 
dial velocity of the satellite. The time difference between 
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FIG. 2. (Color online) Comparison of the measured time of 
detection with the expected time of photon arrival. The com¬ 
patibility of the fit parameters for LEO and MEO satellites 
demonstrates the robustness of the analysis method indepen¬ 
dently from the satellite distance. 


two consecutive Ge/ deviates from the period of 10 ns 
imposed at the generation by a factor 1 -|- 2vnlc, where 
vn is the satellite radial velocity and c is the speed of 
light as shown in Fig. [^bottom). 

In [5D], the Ge/’s were obtained by calculating the 
satellite orbital parameters. Here, similarly to [22]: we 
have introduced an ancillary pulsed beam to simplify the 
prediction of Ge/- Such ancillary beam realize the satel¬ 
lite laser ranging (SLR) in real time and it is derived 
from the same oscillator used for the generation of the 
single photon train, thus providing an intrinsic temporal 
synchronization between the two beams. The SLR beam 
is generated by a three stages laser amplifier system, and 
consists of a train of 100 mJ laser pulses with a repetition 
rate of 10 Hz. It is combined together with the single 
photon beam toward and from the satellite. The exit and 
return instants of the SLR beam are detected using fast 
photodiodes and time-tagged by the same TDG used to 
register the PMT detections. Since both beams follow 
the same path to and from the satellite, the actual tre/ 
of every pulse of the 100 MHz train may be obtained in 
real time and with sub-ns accuracy by the Doppler effect 
measured on the SLR signal (Fig. [^bottom). 

We tested the synchronization method with two satel¬ 
lites in different orbits: Ajisai, a LEO satellite with al¬ 
titude of 1490 km also used in (20] 122] , and LAGEOS- 
2, a MEO satellite with altitude of 5620 km. Eor both 
satellites, we evaluated the difference of the measured 
time-tags of a complete passage with the nearest t^ef 
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calculated by the ancillary SLR pulses. The results are 
presented in Fig. 2 without any selection on the mean 
number of photons hitting the satellite: in both cases, a 
peak exceeding the background was obtained. 

We extracted the residual offset (Ag) and the standard 
deviation {aa) by Gaussian fitting the distribution of the 
differences. The residual offset is Ag = 0.02 ± 0.02 ns in 
both cases and results compatible with 0, thus demon¬ 
strating the validity of this synchronization technique for 
both the LEO and MEO satellites. The measured stan¬ 
dard deviations for the two passages are 0.49 ± 0.02 ns 
for LAGEOS-2 and 0.54 ± 0.02 ns for Ajisai. These two 
values are compatible within the experimental errors, in¬ 
dicating that the synchronization procedure can be ef¬ 
fectively used also for larger satellite distance. In both 
cases, the temporal precision of the single photon counts 
is essentially determined by the 0.50 ns jitter of the PMT 
detectors. 

Link budget analysis - The radar equation [25] predicts 
that the overall link losses, namely the ratio between the 
received {frx) and transmitted {ftx) photon rate, depend 
on the satellite distance R as follows: 



ftx 


aGt 


Tf{R) 

i?4 


( 1 ) 


where a = 237 'wf is a coefficient that includes optical at¬ 
tenuation, geometrical factors, satellite cross section and 
detection efficiency for LAGEOS-2, G* is the transmis¬ 
sion gain and Ta the one-way atmospheric transmission 
coefficient. 

InQ, the detection rate depends on the instantaneous 
satellite slant distance R both directly, due to the factor 
l/i?'*, and implicitly in Ta, which is given by 


2 R (h-t+Rs) 


Ta{R) = 


-(ht+He) 




( 2 ) 


where Tg = 0.89 is the zenith transmission coefficient at 
sea level, /i* = 537 m the observatory altitude above sea 
level. Re = 6371 km the Earth radius, hs = 5620 km the 
satellite altitude and /iscaie = 1.2 km a scaling factor pS] . 
We note that the slant distance R is in general larger than 
the satellite altitude show in Fig. [l] 

We analyzed the entire passage of LAGEOS-2 by di¬ 
viding it in intervals of 60 s each, in which the mean 
return rate (frx) was computed. To point out the events 
to be ascribed to the 100 MHz train from the uncorre¬ 
lated background, we subtracted to the total detection 
rate the background rate measured from the events dis¬ 
tant more than 6 ac form Ge/j where ac is the standard 
deviation of the Gaussian fit of Fig. 

Fig. shows the satellite distance (top) and the frx 
(middle) as a function of the elapsed time from the be¬ 
ginning of the satellite passage. It is possible to no¬ 
tice a drop in frx between minute 17 and 27. This is 
due to a larger telescope pointing error near the Zenith. 
With the exclusion of this interval, frx has been fit 
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FIG. 3. (Color online) Top: satellite distance to the MLRO 
ground station, obtained with time of flight measurements 
of SLR pulses. Middle: single photon detection rate aver¬ 
aged over 60 sec compared with the global fit. The shaded 
area represent the 99.7% of confidence level of the fit. Bot¬ 
tom: detection rate of single photons averaged over 10 sec for 
different values of ^sat: green circles, yellow circles and red 
diamonds respectively correspond to ^sat < 1, 1 < Msat < 2, 
and /isat > 2. 


with Q, leaving as free parameter the transmitter gain 

Gt = jiexp 


-2 ? 


i’t)' 


where 9t is the fixed telescope 


divergence and Op is the time varying pointing error. The 
result of the fit for the whole passage of LAGEOS-2 is 
Gt = (4.7 ± 0.2) X 10®. Excluding the above mentioned 
period, the total photon transmission time was then 32 
min, that greatly extend the link durations previously 
reported pni422j . that were at the level of few tens of 
seconds per satellite passage. This is due to two rea¬ 
sons: the duration of line-of-sight with MEO satellites is 
longer with respect to LEO, and the pointing accuracy is 
increased, due to the lower angular velocity of the satel¬ 
lite, and lower influence of gravitational perturbation of 
the Earth crust. 
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Msat estimation - The mean number of photons hit¬ 
ting the satellite fluctuates due to turbulence and point¬ 
ing errors. This effect is evident in the deviation of the 
measured detection rate from the global fit. To mitigate 
this effect, we based the assessment of Usat not on the 
mean up-link attenuation but as the average number of 
received photon per pulse ^rec divided by the down-link 
part of the radar equation. More specifically ^sat can be 
obtained from 


S 1 

Mrec = t^sat “7 "~Tr ~ a Vdet j (3) 

AcCR P Neff 


where Accr = 11.4 cm^ m is the CCR reflective area, 
E = ISfcm^ the Lageos cross-section, p = 0.89 [3S] the 
CCR reflectivity at normal incidence. At = 1.7357 
the telescope area, rjrx = 0.1306 the transmission of the 
whole receiving apparatus, rj^et =0.1 the detector ef¬ 
ficiency, and Neff = 9.88 [37] the effective number of 
CCR, averaged over all orientation. 

To reduce the effect of the fluctuation of the point¬ 
ing error on p^at estimation, we divided the whole pas¬ 
sage in intervals of 10 s, calculating for each interval 
the mean prec- With this method we computed psat for 
each analysis slice, distinguishing those where psat < 1, 

1 < Psat < 2, and psat > 2 (Fig. bottom). The first 
two cases have a p that can be considered practical to be 
used in protocols such as BB84 with decoy states [331135] . 
As shown in Fig.[^ we received pulses with less than one 
photon per pulse when the satellite is at the largest dis¬ 
tance from the ground station, corresponding to the max¬ 
imum attenuation, and when the telescope had pointing 
errors, which increases the up-link losses. 

By considering the intervals corresponding to p^at ^ 1 
it has been possible to obtain a total integration time 
of 200 s. For these periods we compared the times¬ 
tamp of each photon detection with the nearest tref, with 
the same technique devised for the analysis reported in 
Fig.H As shown in the upper graph of Fig. [^ we ob¬ 
tained a peak with statistical significance of (5.1 ± 1.2 )ct 
over background (here a is the standard deviation of the 
background counts). The mean signal intensity for this 
intervals is psat = 0.55 ± 0.06 while the mean down¬ 
link attenuation is 72.3 dB. The average detection rate is 
3.0 ± 0.3 counts js, and the peak signal to noise ratio is 
SNR= 1.5 ±0.1. 

We repeated the same analysis for intervals with psat ± 

2 (Fig. ID bottom). This resulted in an integration pe¬ 
riod of 510 s, a statistical significance of the peak over 
background of (15.1 ± 1.3 )(t, a mean psat = 1-19 ± 0.05 
and a similar downlink attenuation, of 72.4 dB. The 
mean detection rate is 6.8 ± 0.3 counts/s and the peak 
SNR = 1.8 ±0.1. 

Based of these observations, we extrapolate the feasi¬ 
bility of QC with MEO satellites by simulating the pho¬ 
ton link budget using the best current technology of de¬ 
tectors (silicon based single photon avalanche detectors) 




FIG. 4. (Color online) Histograms of detections as a function 
of deviation A from expected time pet for different values of 
Psat- When only intervals with paat < 1 are selected (top), 
a peak is observed at A = 0 with statistical significance over 
background of 5.1 cr. If the analysis is extended to paat < 
2 (bottom) the peak significance reaches Bin size is 

0.4 ns. Robustness of the peaks has been verified by changing 
the bin size. 


and under the hypothesis of effectively coupling on them 
the satellite returns. From Ref. [30], peak detection effi¬ 
ciency is taken as ps = 48%, FWFIM timing accuracy as 
50 ps and dark count rate as DCs — 350 counts/s. These 
values significantly overpass the ones available for present 
experiment, ppmt = 10%, FWHMpmt = 1-22 ns and 
DCpMT = 50 counts/s, that are relative to the PMT 
technology, which was imposed by the current optical 
system coupled to the Coude path. We note that 50 ps 
are compatible with the expected temporal broadening 
of the pulse due to propagation, that are of the order of 
few ps ED. Moreover, by adopting a modulated CCR 
as proposed in [33] with a light amplitude modulator al¬ 
lowing for pfixed = 0.6 during the whole satellite pas¬ 
sage, it would be possible to enhance the detection rate 
by a factor = 5.2 while reducing the dark 

count rate within the coincidence window by a factor 

PCpMT FWHMpmt ^ o At. 

DCs FWHMs — 

Because of these steps, the SNR increases by a factor 
18 with respect to the values here obtained, correspond¬ 
ing to a quantum bit error ratio QBER ~ 3.6%. In the 
hypothesis that no other effects would step in the bud¬ 
get [53], QKD from a MEO satellite appears as feasible. 
For extending QC to higher orbit such as the GNSS, be¬ 
tween 19000 km and 23000 km of altitude, a possible 
improvement of the detector dark counts down to ~ 100 
counts/s is necessary. In this case the signal would be 
reduced by a factor 10.8, due to the 1/i?^ factor in (§, 
giving a value of QBER ~ 6.6% that is compatible with 
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the implementation of QKD for such GNSS satellites. 

A possible improvement could come from the use of 
superconducting single photon detectors. With nanowire 
superconducting technology it is possible to obtain effi¬ 
ciency up to 80%, dark count rates as low as 10 Hz and 
time jitter of 40 ps. With these parameters it would 
be possible to obtain a QBER = 0.5% form a satel¬ 
lite with the same orbital parameters as LAGEOS, and 
QBER = 3% from a GNSS satellite. However it mus be 
noticed that the smaller sensitive area of superconductive 
devices, which reaches ~ 16^m compared to ~ 200/iTO of 
silicon based detectors, might lead to a lower efficiency 
due to the optical coupling. 

Conclusions - We report on the experimental exchange 
of single photons from the LAGEOS-2 satellite to the 
MLRO ground station, with a link length exceeding 
7000 km. The photon transmission has been kept sta¬ 
ble for most of the passage of the MEO satellite, allow¬ 
ing for the assessment of the signal attenuation due to 
beam divergence and telescope pointing error by means 
of the established link models. By using an active source 
on satellite and commercially available silicon detectors, 
a QBER suitable for the implementation of a QKD link 
from a MEO satellite appears as feasible. This also paves 
the way to endeavoring QG with GEO satellites and even¬ 
tually to fundamental tests of combined Quantum Me¬ 
chanics and General Relativity. 
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